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ABSTRACT 

An experimental search for variation in the fundamental coupling constants is strongly 
motivated by modern high-energy physics theories. Comparison of quasar absorption 
s ! ' line spectra with laboratory spectra provides a sensitive probe for variability of the fine 

structure constant, a, over cosmological time-scales. We have previously developed and 
■ applied a new method providing an order of magnitude gain in precision over previous 

optical astrophysical constraints. Here we extend that work by including new quasar 
spectra of damped Lyman-a absorption systems. We also re-analyse our previous 
O !' lower redshift data and confirm our initial results. The constraints on a come from 

simultaneous fitting of absorption lines of subsets of the following species: Mgi, Mgn, 
Alii, Aim, Sin, Cm, Fell, Mil and Znii. We present a detailed description of our 
methods and results based on an analysis of 49 quasar absorption systems (towards 
28 QSOs) covering the redshift range 0.5 < z < 3.5. There is statistical evidence for 
a smaller a at earlier epochs: Aa/a = (—0.72 ± 0.18) x 10~ 5 . The new and original 
samples are independent but separately yield consistent and significant non-zero values 
of Aa/a. We summarise the results of a thorough investigation of systematic effects 
5.^ ■ published in a companion paper. The value we quote above is the raw value, not 

corrected for any of these systematic effects. The only significant systematic effects so 
far identified, if removed from our data, would lead to a more significant deviation of 
Aa/a from zero. 

Key words: atomic data - line: profiles - methods: laboratory - techniques: spec- 
troscopic - quasars: absorption lines - ultraviolet: general 
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1 INTRODUCTION 

High resolution spectroscopy of QSO absorption systems 
provides a precise probe of fundamental physics. In pre- 
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vious papers (Dzuba, Flambaum & Webb 1999a,b; Webb 
et al. 1999, hereafter W99) we introduced and applied a 

new and highly sensitive method for constraining space- 

2 

time variations of the fine structure constant, a = f-. Us- 

' he 

ing 30 Keck/HIRES absorption spectra containing Mg n and 
Fell, we reported a tentative detection of a variation in a: 
Aa/a = (-1.09 ± 0.36) x 10~ 5 , where Aa/a is denned as 

Aa/a = (a, - a )/a , (1) 

ao is the present day value of a and a z is the value at the 
absorption redshift, z. The absorption redshifts ranged from 
~0.6-1.6. 
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Here we continue this work, presenting new measure- 
ments of Aa/a in high redshift damped Lyman-a systems 
(DLAs) over the range 0.9 < z < 3.5, and a re-analysis of 
the W99 data set. Four papers present our new work on this 
subject. A brief summary of all results is given in Webb et 
al. (2001). The present paper provides a detailed description 
of our methods and results. A companion paper (Murphy et 
al. 2001b, hereafter MOlb) presents a thorough examination 
of systematic effects. An additional paper (Murphy et al. 
2001c) presents an analysis of a sample of Siiv QSO ab- 
sorbers. 

The present paper is organized as follows. Section ^ 
describes theoretical motivations for a search for varying 
a. In Section ^ we summarise previous experimental con- 
straints on varying a. Section ^| presents the theoretical ba- 
sis of the new method for determining Aa/a using QSO 
absorbers. Section |^ describes the experimental components 
of our study: high-resolution QSO observations and precise 
laboratory rest wavelength measurements. In Section |^ we 
discuss our analysis technique in detail and present our re- 
sults. For completeness we also summarise the systematic 
effects discussed in detail in MOlb. Section summarises 
our new results and outlines future work. 



2 THEORETICAL BACKGROUND 
2.1 Early historical interest 

The constancy of the fundamental constants began to be 
questioned when Milne (1935, 1937) and Dirac (1937) sug- 
gested that the Newtonian gravitational constant, G, may 
vary in cosmological time. Milne proposed a gravitational 
theory in which different physical 'clocks' - 'ticking' ac- 
cording to different physical processes - 'tick' at differ- 
ent rates as time progresses. The possible consequences of 
this hypothesis for biological and geological history were 
seized upon by Haldane (1937a,b) to explain aspects of 
their discontinuous development, whereas Kothari (1938) 
and Chandrasekhar (1939) pointed out some astronomical 
consequences of Dirac's 'large numbers hypothesis' which 
might permit it to be tested by observation. Dirac's Large 
Numbers Hypothesis attracted considerable interest (e.g. 
Teller 1948; Dicke 1957, 1961; Barrow 1990; Barrow & Tipler 
1996) and the desire to place his idea of varying G on a firmer 
theoretical footing led to the development of scalar-tensor 
generalisations of Einstein's theory of general relativity in 
which the variation of G was described self-consistently by 
the propagation of a scalar field which also acted as a source 
of space-time curvature (Brans & Dicke 1961). 

Jordan (1937, 1939) was the first to consider how 
Dirac's hypothesis might be applied to forces other than 
gravity, but rejected the possibility of time variations in the 
weak interaction or the electron-proton mass ratio (Jordan 
1937, 1939). Later, Gamow (1967) showed that geological 
problems caused by Dirac's proposed variation in G that 
had been pointed out by Teller (1948) could be avoided 
by interpreting the large numbers hypothesis as requiring 
the time-variation of the charge on the electron, e, rather 
than G. Observational limits on possible variations of a or 
the weak and strong interactions were more restrictive than 
those on G although no self-consistent theory permitting 



the variation of non-gravitational force constants was de- 
veloped to make the observational limits rigorous. Landau 
(1955) proposed that variation in a could be connected to 
renormalization rules in quantum electrodynamics. 

Current motivation for a search for a varying a comes 
from modern theories which attempt to unify gravity with 
other interactions. 



2.2 Modern theoretical motivations for varying a 

At present there are many theoretical motivations for a 
search for varying a. However, these theories, being rel- 
atively recent, may not be widely appreciated within the 
context of astronomical tests of their viability. We therefore 
summarise recent ideas below. 

2.2.1 Multi- dimensional unification theories 

One of the greatest challenges facing modern theoretical 
physics is the quantization of the gravitational interaction. 
Early attempts at creating such a unification were made us- 
ing a geometrization in (4 + D)-dimensional curved space- 
times in the spirit of the old Kaluza-Klein scenario to unite 
gravity and electromagnetism (Kaluza 1921; Klein 1926). 
Three-dimensional gauge couplings like a, a wca k, and a a tron g 
then must vary as the inverse square of the mean scale of the 
extra dimensions. Thus, the evolution of the scale size of the 
extra dimensions is related to variability of the low-energy 
coupling constants in the 4-dimensional subspace in simple 
Kaluza-Klein and superstring theories. 

Damour & Polyakov (1994) showed that cosmological 
variation in a may proceed at different rates at different 
points in space-time (see also Forgacs & Horvath 1979; 
Barrow 1987; Li & Gott 1998). Various functional forms 
for time variation of a/G have been derived using Kaluza- 
Klein theory and the assumption of constant masses. Cho- 
dos & Detweiler (1980) find a/G oc t and Freund (1982) 
finds a/G oc t 1/A . Wu & Wang (1986) have found that 
oc H ~ — 10 _11±1 yr _1 where the proportionality constant 
is unknown. It is interesting to note that using a similar 
analysis to Wu & Wang, Maeda (1988) derived a oc i" 4/3 . 

Marciano (1984) discussed the self-consistency relations 
required if there are simultaneous variations of different con- 
stants in unified gauge theories. He also examined any pos- 
sible non-monotonic variation in a with t, using the running 
coupling dependence of strong, weak, and electromagnetic 
interactions to produce self-consistent predictions for the si- 
multaneous variation of more than one coupling or mass 
ratio. These were discussed in Drinkwater et al. (1998). Typ- 
ically, variation in G and a could be linked by relations of 
the form ^ ~ 

These features of extra dimensions mean that in all 
string theories (and M-theory, of which they are presumed to 
be limiting cases, Hofava & Witten 1996a, b), any extra di- 
mensions of space need to be held static with high accuracy 
in order to avoid gross conflict with observation. In the cur- 
rently popular scenarios for M-theory (e.g. Antoniadis et al. 
1998; Arkani-Hamed, Dimopoulos & Dvali 1998; Randall & 
Sundrum 1999a,b), only the gravitational force is assumed 
to act in all (> 3) spatial dimensions (the 'bulk') whilst 
all other interactions act only in 3-dimensional space (the 
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'brane'). Observations of the constancy of 3-dimensional 
non-gravitational constants governing the strengths of those 
interactions in three dimensions (like a) could therefore be of 
crucial importance in testing this theoretical scenario. The 
scale of the extra dimensions was until recently believed to 
be necessarily of order the Planck scale ~ 10 -33 cm but it 
is possible they could be as large as ~ 0.01 mm leading to 
large changes in the form of the Newtonian law of gravity 
on length scales below this distance. At present, this is con- 
sistent with all known experiments. The challenge of testing 
the law of gravity on sub-millimetre scales is considerable 
because other forces strongly dominate the effects of gravity 
on these scales. 



2.2.2 Scalar Theories 

The first self consistent theory of electromagnetism which in- 
corporates varying a, and which reduces to Maxwell's theory 
in the limit of constant a, was that developed by Bekenstein 
(1982). It introduces a scalar field whose variation in space 
and time produces effective variation in the electron charge 
(or, equivalently, the permittivity of free space). The form of 
the propagation equation for the evolution of the scalar field 
is strongly constrained by the natural requirements that it 
be second order, causal and linearly coupled to the matter 
density. This particular theory has been constrained by sev- 
eral astronomical observations (Livio & Stiavelli 1998) and 
used to illustrate consequences of space and time variations 
of constants (Barrow & O'Toole 2000). 

The general structure of Bekenstein's theory can also 
be applied to other couplings, and it has recently been used 
to investigate the problem of varying the strong coupling 
constant of QCD by Chamoun, Landau & Vucetich (2000). 
Bekenstein makes the simplifying assumption that the grav- 
itational sector of the theory be unchanged from that of 
general relativity. This amounts to neglecting the contri- 
bution of the kinetic energy of the scalar field associated 
with a variation to the expansion dynamics of the universe. 
This defect can be avoided in other versions of this the- 
ory (Barrow & Magueijo 2000; Magueijo 2000). It should 
be emphasised that there is no single self-consistent theory 
incorporating varying a and theoretical limits need to be 
quoted in conjunction with the theoretical framework used 
to confront observational data. If different interactions are 
unified by electro-weak or grand unified theories then the 
ramifications of introducing varying a will be deeper and 
wider and a scalar theory would need to accommodate the 
effects of varying a on the other gauge couplings. 

2.2.3 Varying speed of light theories 

One appeal of varying speed of light theories is that they 
provide alternative potential solutions to a range of cosmo- 
logical problems such as the flatness, horizon and monopole 
problems (Moffat 1993; Barrow & Magueijo 1999a,b; Al- 
brecht & Magueijo 1999; Barrow 1999). Barrow & Magueijo 
(1998) have also expressed Bekenstein's (1982) varying-e 
theory in terms of varying c. 

Barrow & Magueijo (2000) developed a particular the- 
ory for varying c (or a) in which the stress contributed by 
the cosmological constant varies through the combination 



Ac 2 . They also showed how the observed non-zero cosmo- 
logical acceleration (Schmidt et al. 1998; Perlmutter et al. 
1999) might be linked to a varying a. They examined a class 
of varying-c theories in which changes in c are driven by a 
scalar field which is coupled to the gravitational effect of 
pressure. Changes in c convert the A energy density into ra- 
diation, preventing A from dominating the universe in the 
radiation epoch. As pressureless matter begins to dominate, 
variations in c slow, and A begins to dominate. The very slow 
variation of the scalar field driving changes in c is like that 
of scalar fields in quintessence theories in the presence of 
perfect fluids. This type of theory allows the time variations 
in c or a to be ~ 10 -5 Ho at z ~ 1 and yet the associated A 
term can be dominant today and produce acceleration (the 
10~ 5 factor is created naturally in this theory by the ratio 
of the radiation to matter density in the universe today). 
We therefore expect only very small residual variation in a 
in a A dominated universe. The existence of simple theories 
of this sort provides a strong motivation for studies like the 
present one. 

2.2.4 Other theories 

Several authors have suggested other motivations for vary- 
ing a. For example, Hill & Ross (1988a,b) (see also Hill, 
Steinhardt & Turner 1990) have suggested that oscillatory 
variations in a might be expected if there exists a soft bo- 
son field that has electromagnetic couplings to either the 
Maxwell scalar, F^F 11 " , or to the electron state, ^ e ipe. If 
the boson field has a mass m then the frequency of emis- 
sion and absorption of radiation will vary with a frequency 
~ m~ . An oscillatory a might also be attributable to the 
decaying oscillations of a scalar field that couples to F^F^" . 
Such a theory is allowed by enforcing an approximate global 
symmetry (Carroll 1998). 

In summary, there is much current theoretical motiva- 
tion for experimental searches for variability in a, particu- 
larly at medium to high redshifts. A detection of variation 
may allow us to probe the extra dimensions expected from 
many modern unified theories. A search for varying a pro- 
vides an important tool for seeking out and constraining the 
form of new physics. 



3 PREVIOUS EXPERIMENTAL WORK 

3.1 Terrestrial and laboratory constraints 

Direct laboratory measurements provide interesting con- 
straints on a time-varying a. Prestage, Robert & Maleki 
(1995) introduced a new technique comparing the rates of 
clocks based on ultra-stable oscillators. Relativistic correc- 
tions are order (Z„a) 2 where Z„ is the atomic number or 
nuclear charge. By comparing the rates of two clocks based 
on different atoms (H-maser and Hg + ) over a 140 day pe- 
riod, they were able to constrain \a/a\ < 3.7 x 10~ 14 yr -1 
(i.e. \Aa/a\ < 1.4 x 10~ 14 ). 

Model dependent upper limits on any variation have 
been claimed from an analysis of the Oklo phenomenon - a 
natural nuclear fission reactor that operated in Gabon, West 
Africa, ~1.8 billion years ago. Shlyakhter (1976) originally 
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constrained the energy of the nuclear resonance level in the 
150 Sm isotope and related this to upper bounds on \Aa/a\. 
Damour & Dyson (1996) have also analysed the problem 
with different model assumptions and claim —0.9 x 1(F 7 < 
Aa/a < 1.2 x l(T r . Recently, Fujii et al. (2000) obtained 
somewhat tighter constraints using new samples from the 
Oklo reactor: Aa/a = (-0.04 ±0.15) x 10~ 7 . 

However, the Oklo limit corresponds to variations at 
very low 'redshift', z ~ 0.1, in a non-cosmological environ- 
ment. In the absence of a detailed theory giving both time 
and space variations of a, it is dangerous to compare time 
variations in areas of different gravitational potential (Bar- 
row & O'Toole 2000). Also, there are considerable uncer- 
tainties regarding the complicated Oklo environment and 
any limits on a derived from it can be weakened by allowing 
other interaction strengths and mass ratios to vary in time 
as well. 

3.2 Constraints from QSO spectra 

Savedoff (1956) first analysed doublet separations seen in 
galaxy emission spectra to constrain possible time variation 
of a. Absorption lines in intervening clouds along the line 
of sight to QSOs are substantially narrower than intrinsic 
emission lines and therefore provide a more precise probe 
of a. This prompted Bahcall, Sargent & Schmidt (1967) to 
use alkali-doublet spacings seen in the absorption spectra of 
QSOs to arrive at firm upper limits on any variation. The 
appeal of this method is the huge time span probed by high- 
redshift observations. Wolfe, Brown & Roberts (1976) first 
applied the alkali-doublet (AD) method to intervening Mg II 
absorption. Varshalovich, Panchuk & Ivanchik (1996) have 
used the AD method for many absorbers to obtain upper 
limits on any variation at redshifts z ~ 2.8-3.1: Aa/a — 
(2± 7) x 10~ 5 . 

Varshalovich, Potekhin & Ivanchik (2000) have recently 
analysed 16 Siiv absorption systems (towards 6 QSOs) us- 
ing the AD method. They find a mean Aa/a — (—4.6 ± 

4.3 ± 1.4) x 10~ 5 using the Siiv laboratory wavelengths 
of Morton (1991). The systematic component to the error 
above (1.4 x 10 ) allows for the relative uncertainties of the 
Si IV laboratory wavelengths. However, Griesmann & Kling 
(2000) have increased the absolute precision of these wave- 
lengths and their new values indicate that the Varshalovich 
et al. (2000) result should be corrected by « +11.3 x 10" 5 to 
Aa/a ~ (+7 ±4) x 10~ J suggesting that the systematic er- 
ror component added by Varshalovich et al. was significantly 
underestimated. 

In Murphy et al. (2001c) we have used the new Gries- 
mann & Kling laboratory wavelengths to constrain Aa/a 
using 21 Siiv doublets towards 8 QSOs. The spectral res- 
olution and SNR are higher than for the Varshalovich et 
al. sample. This, and the higher precision laboratory wave- 
lengths, provide a more stringent constraint of Aa/a — 
(—0.5 ± 1.3) x 10~ 5 . This result is nevertheless a factor of 7 
less precise than the results we present in this paper, high- 
lighting the strength of the methods we use. 

The optical techniques above were based on the AD 
method: constraining a from the splitting of single ADs such 
as Siiv, Civ and Mgn. Very recently, however, Dzuba et 
al. (1999a) and W99 proposed a new method that allows 
an order of magnitude increase in precision. The method 



is based on measuring the wavelength separation between 
the resonance transitions of different ionic species with no 
restriction on the multiplet to which the transitions belong. 
We outline this many-multiplet (MM) method and discuss 
its advantages in Section ^| of this paper. 

Radio spectra of QSO absorption clouds also provide 
the opportunity for tightly constraining possible variations 
in the fundamental constants. Wolfe et al. (1976) demon- 
strated that combining optical resonance and Hi 21 cm 
hyperfine absorption spectra could yield constraints on 
a 2 g p ^f-. Here, g p is the proton g- factor and m e and m p 
are the mass of the electron and proton respectively; see 
also the discussion by Pagel (1983). Varshalovich & Potekhin 
(1996) used a comparison of H 1 21 cm hyperfine and molec- 
ular rotational absorption to constrain variations in m p but 
Drinkwater et al. (1998) showed that such a combination ac- 
tually constrains y = a 2 g p . In Murphy et al. (2001d) we have 
also made a more complete analysis of the two Drinkwater 
et al. radio spectra to obtain constraints on y. Assuming g p 
to be constant, we find Aa/a = (-0.10 ± 0.22) x 10~ 5 and 
Aa/a = (-0.08 ± 0.27) x 10~ 5 at z = 0.25 and 0.68 respec- 
tively. If we note the low-z points in the lower panel (binned 
data) of Fig. 3 then we see that our results are also consis- 
tent with the two radio points. Carilli et al. (2000) have also 
recently analysed independent spectra of these two QSOs 
and find an upper limit on any variation consistent with our 
measurements: |Aa/a| < 0.85 x 10~ 5 (again assuming g p to 
be constant). 

Cowie & Songaila (1995) have also obtained constraints 
at z ~ 1.78. They compare the observed frequencies of the 
H 1 21 cm and UV C I absorption lines which are sensitive to 
variations in g p a 2 ^-. If one assumes that g p — is constant 
then their constraint corresponds to Aa/a = (0.35 ±0.55) x 
10~ 5 . 



3.3 Early universe (CMB and BBN) 

A smaller value of a in the past would change the ioniza- 
tion history of the universe, postponing the recombination 
of electrons and protons, i.e. last-scattering would occur at 
lower redshift. It would also alter the ratio of baryons to pho- 
tons at last-scattering, leading to changes in both the am- 
plitudes and positions of features in the cosmic microwave 
background (CMB) power spectrum, primarily at angular 
scales Si 1° (Hannestad 1999; Kaplinghat, Scherrer & Turner 
1999; Avelino, Martins & Rocha 2000; Battye, Crittenden & 
Weller 2000; Kujat & Scherrer 2000). 

Kolb, Perry & Walker (1986) have discussed how limits 
can be placed on Aa/a at the time of big bang nucleosyn- 
thesis (BBN) (corresponding to a redshift z ~ 10 8 — 10 9 ) 
by assuming a simple scaling between the value of a and 
the neutron-proton mass difference. Using observations of 
the 4 He abundance they derive limits of \a/a\ < 15 x 
10~ 15 /ii 00 yr~\ corresponding to \Aa/a\ < 9.9 x 10~ 5 . 

Barrow (1987) and Campbell & Olive (1995) consid- 
ered the simultaneous variation of weak, strong and electro- 
magnetic couplings on primordial nucleosynthesis. However, 
these limits are beset by a crucial uncertainty as to the elec- 
tromagnetic contributions to the neutron-proton mass dif- 
ference. Much weaker limits are possible if attention is re- 
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stricted to the nuclear interaction effects on the nucleosyn- 
thesis of deuterium, helium-3, and lithium-7. 

Bergstrom, Iguri & Rubinstein (1999) have re-examined 
all light elements up to 7 Li. By including elements heavier 
than 4 He they avoid the problem of the unknown electro- 
magnetic contributions to the neutron-proton mass differ- 
ence, and argue that a more realistic limit, given the present 
observational uncertainties in the light element abundances 
is more like 2 x 10 -2 . 



4 THEORETICAL BASIS OF THE 
MANY-MULTIPLET METHOD 

The MM method allows an order of magnitude greater pre- 
cision compared to the previous AD method. The full details 
behind this technique are presented in detail in Dzuba et al. 
(1999a,b, 2001) and W99 and so we present only the salient 
features here. 

To explain the MM method, we begin with a simple 
analytic approach to estimate the relativistic effects in tran- 
sition frequencies. If we consider a many-electron atom/ion 
then the relativistic correction, A, to the energy of the ex- 
ternal electron can be written as 



A oc (Z n af \E 



3/2 



3 + 1/2 



C(j,l) 



(2) 



where Z n is the nuclear charge, E is the electron energy 
(E < 0, \E\ is the ionization potential) and j and I are the 
total and orbital electron angular momenta. The contribu- 
tion to the relativistic correction from many-body effects is 
described by C(j, I). For s and p orbitals, C(j, I) « 0.6 and 
is of similar magnitude for d orbitals. Equation [| therefore 
provides a general strategy for probing the relativistic cor- 
rections in resonance transitions. 

For example, consider comparison of the transition en- 
ergies of two s-p transitions, one in a light ion, the other 
in a heavy ion (i.e. low and high Z n respectively). The Z„ 
term dominates so the relativistic corrections to the tran- 
sition energies will differ greatly. Thus, comparison of the 
spectra of light and heavy species is a sensitive probe of a. 

As a further example, consider an s—p and a d-p transi- 
tion in a heavy species. The corrections will be large in each 
case but will be of opposite sign since the many-body cor- 
rections, C(j,l), begin to dominate with increasing j. This 
situation also allows tight constraints to be placed on a. 

Thus, comparing spectra of transitions from different 
multiplets and different atoms or ions, provides a sensitive 
method for probing variations in a. In comparison, the fine 
splitting of an s-p doublet will be substantially smaller than 
the absolute shift in the s-p transition energy since the ex- 
cited p electron, with relatively small \E\, will have much 
smaller relativistic corrections than the s electron. There- 
fore, the AD method is relatively insensitive to variations in 
a. 

More formally, the energy equation for a transition from 
the ground state, within a particular multiplet, at a redshift 
z can be written as 



E z — E c + QiZ n 



{-)' 
\a J 



+ 



K^Zlf^) +K 2 (LS) 2 Z 4 J^) , 
V Qo / \ ao J 

© 2000 RAS, MNRAS 000, §-0 



(3) 



where a z may or may not be equal to the laboratory value, 
ao. Here, L and S are the electron total orbital angular mo- 
mentum and total spin respectively and E c is the energy of 
the configuration centre. Qi, Ki and K 2 are relativistic co- 
efficients which have been accurately computed in Dzuba et 
al. (1999a,b, 2001) using ab initio many-body calculations. 
These calculations were tested by comparison with exper- 
imental data for energy levels and fine structure intervals. 
Equation | forms the basis of the MM method since it can 
be applied to any resonance transition observed in the QSO 
spectra - another advantage it has over the AD method. 

For our purposes, the most convenient form of equation 
|§] is written as 



to z — uj + qix + q 2 y ■ 



(4) 



where to z is the wavenumber in the rest-frame of the cloud, 
at redshift 2, Do is the wavenumber as measured on Earth 
and x and y contain the information about a possible non- 
zero Act /a: 



(^) 2 -l and „ S (2-) 4 - 1 
V ao / \ ao J 



(5) 



The values of the relativistic corrections, qi and 52, are pre- 
sented in Table [j] for all transitions commonly seen in opti- 
cal QSO spectra. The values of the qi coefficients are typi- 
cally an order of magnitude larger than the q 2 coefficients. 
Therefore, it is the relative magnitudes of qi for different 
transitions that characterize our ability to constrain Aq/q. 
The uncertainty in the qi co-efficients varies from ion-to-ion 
and is ~ 5 cm -1 for alkali-like ions such as Mgn. However, 
the co-efficients for the compilcated Nill transitions carry 
a larger uncertainty ~ 200 cm -1 (Dzuba et al. 2001). The 
form of equation ^] is very convenient since the second and 
third terms only contribute if Aa/a 7^ 0: errors in the qi and 
52 coefficients cannot lead to an artificial non-zero value of 
Aa/a. 

From equation |^ we can see that qi will be much larger 
in heavier ions due to the Z\ dependence. Indeed, one can 
see in Table [j] several combinations of transitions from dif- 
ferent ions, of very different mass, showing large differences 
between their respective qi coefficients. The combination of 
the Mg 1 A2853 line, the Mg 11 doublet and the five Fe 11 lines 
between A2344 and A2600 was used in W99 since the Fell 
qi coefficients typically differ from those of the Mg lines by 
almost an order of magnitude. It is this comparison of spec- 
tra from different ions that allows the order of magnitude 
increase in precision over previous methods. For compari- 
son, we present the line shifts in wavenumber, wavelength 
and velocity space arising from a shift of Aa/a = +10~ 5 in 
Table |. 

Further examination of the qi coefficients in Table [l] 
reveals several other very useful combinations of lines. In 
particular, note the large negative values for the Cm lines. 
A comparison between Cr 11 spectra and that of Zn 11 - the 
A2026 transition having the largest magnitude value of qi - 
provides the most sensitive combinations for probing non- 
zero Aa/a. Also, note the coefficients for Nin. Here, within 
the same species (but different multiplets), we have both 
positive and negative values for qi . This unique case is due to 
the very complicated multiplet structure of Ni 11 as discussed 
in detail in Dzuba et al. (2001). 
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Table 1. Atomic data for those transitions used in our analysis. The origins of the composite wavenumbers (u>o) and wavelengths (Ao) 
are detailed in Section 5.2. We give the ground and excited state configurations in columns 4 and 5. The ID letters in column 6 are used 
in Table Q to indicate the transitions used in our Voigt profile fits to each absorption system. The seventh column shows the ionization 
potential (IP) a and the eighth shows the literature oscillator strengths, /. Values of qi and q 2 are also shown. The Sin, Alii and Aim 
wav enumbers have been scaled from their literature values due to the Norlen/Whaling et al. calibration difference described in Section 
5.2. The isotopic structures for the Mg and Si transitions are given in Table 



Ion 




u 


Ground 


Upper 


ID 


IP 


f 


qi 


92 




A 


cm -1 


state 


state 




cV 




cm -1 


cm -1 


Ms 1 


2852.96310(8) 


35051. 277(1) 6 


3s 2 % 


3s3p 1 Pi 


■ c \ 




1.810 9 


106 


-10 


Mg 11 


2796 S54ff2l 


35760 f,Af,(1\ b 


3s 2 Si /o 


3n 2 Po In 


b 


7.7 


0.6295' 1 


211 







2803.5315(2) 


35669. 298(2) 6 




3r> 2 Pi m 
* J i J L 1/2 


Q 




0.3083' 1 


120 





AID 


1670.7887(1) 


59851.972(4) c 


3s 2 ^0 


3s3p 1 Pi 


d 


6.0 


1.88 9 


270 





Aim 


1854.71841(3) 


53916.540(l) c 


3s 2 S 1/2 


3p 2 P 3 /2 


c 


18.9 


0.2689 


464 







1862.79126(7) 


53682.880(2) c 


3p 2 Pi/ 2 


f 




0.5399 


216 





Sin 


1526.70709(2) 


65500.4492(7) c 


3s 2 3p 2 P° 


3s 2 4s 2 S 1/2 


g 


8.2 


0.116' 


24 


22 




1808.01301(1) 


55309.3365(4) c 


3s3p 2 2 D 3/2 


h 




0.00218 9 


525 


3 


Cr 11 


2056.25693(8) 


48632.055(2) d 


3d 5 6 S 5/2 


3d 4 4p 6 P 7/2 


i 


6.8 


0.105 J 


-1030 


-13 




2062.23610(8) 


48491.053(2) d 


3dHp 6 P 5/2 


j 




0.078 J 


-1168 


-16 




2066.16403(8) 


48398.868(2) d 




3dHp 6 P 3/2 


k 




0.0515 J 


-1267 


-9 


Fen 


1608.45085(8) 


62171. 625(3) e 


3d 6 4s a 6 D 9/2 


3d 6 4p y 6 P r/2 


1 


7.9 


0.0619 fe 


1002 


141 




1611.20034(8) 


62065.528(3) e 


3dHp y 4 F 7/2 


m 




0.00102 fc 


1110 


48 




2344.2130(1) 


42658.2404(2)-^ 




3d 6 4p z 6 P 7/2 


11 




0.110 fc 


1325 


47 




2374.4603(1) 


42114.8329(2)-^ 




3d 6 4p z 6 F g/2 







0.0326 fe 


1730 


26 




2382.7642(1) 


41968.0642(2)-^ 




3d 6 4pz 6 F 11/2 


P 




0.300* 


1580 


29 




2586.6496(1) 


38660.0494(2)-'' 




3d 6 4p z 6 D° /2 


q 




0.0684 fe 


1687 


-36 




2600.1725(1) 


38458.9871(2)-/' 




3d 6 4p z 6 D° /2 


r 




0.213 fc 


1449 


2 


Ni 11 


1709.6042(1) 


58493.071(4) d 


3d 9 2 D 5/2 


3d 8 4pz 2 F 5/2 


s 


7.6 


0.0348' 


800 







1741.5531(1) 


57420.013(4) d 


3d 8 4p 2 2 D 5/2 


t 




0.0419' 


-700 







1751.9157(1) 


57080.373(4) d 




3d 8 4p z 2 F 7/2 


u 




0.0264' 


-300 





Zn II 


2026.13709(8) 


49355.002(2) d 


3d 10 4s 2 S 1/2 


3d 10 4p 2 P 3/2 
3d 10 4p 2 P 3/2 


V 


9.4 


0.489 3 


2291 


94 




2062.66045(9) 


48481.077(2) d 


w 




0.256-? 


1445 


66 



a IP is defined here to be the energy required to form the ion in question from the ion with a unit lower charge; 'Pickering, Thorne & 
Webb (1998); c Griesmann & Kling (2000); d Pickering et al. (2000); e S. Johansson (private communication); ^Nave et al. (1991); 9Morton 
(1991); ^Verner, Verner & Ferland (1996); ^Dufton et al. (1983); ^Bergeson & Lawler (1993); fc Cardelli & Savage (1995); 'Fedchak & 
Lawler (1999) 



We may summarise the advantages of the MM method 
over the AD method as follows: 

• By including all relativistic corrections (i.e. including 
those for the ground state) there is a sensitivity gain of 
around an order of magnitude compared to the AD method. 

• In principle, all transitions appearing in a QSO absorp- 
tion system may be used. This provides an obvious statistical 
gain and a more precise constraint on Aa/a compared to 
using a single AD alone. 

• A further advantage of using many transitions is that 
the velocity structure is determined with much greater reli- 
ability due to the larger range of line strengths (see Section 



5.1 



for further explanation). 
• A very important advantage is that comparison of tran- 
sitions with positive and nega tive q i coefficients minimizes 
systematic effects (see Section 5.1.3 in particular). 



Finally, we note that the q\ and q 2 coefficients for 
A1IIA1670 and SinA1526 in Table |l| have not been pre- 
sented elsewhere before. In the case of AlllA1670, the co- 
efficients were calculated using the same methods described 
in Dzuba et al. (1999a,b). The coefficients for SillA1808 do 
appear in Dzuba et al. (1999b) but the calculation has been 
re-performed to include other Sill lines, including A1526. 



Thus, there is a small (« 5%) difference between the previ- 
ous value and the one presented in Table [j]. This small dif- 
ference is within the precision of the Dzuba et al. (1999a,b) 
calculations and does not significantly affect our determina- 
tion of Aa/a. 



5 DATA REQUIREMENTS, ACQUISITION 
AND REDUCTION 

The fitting of high-resolution (FWHM - 7kms-\ R ~ 
45000), moderate signal-to- noise ratio (SNR ~ 30 per pixel) 
QSO absorption lines routinely yields redshifts of single ve- 
locity components measured to precisions of a(z) ~ 3x 10~ 6 
or o(v) ~ lkms -1 in velocity space (Outram, Chaffee & 
Carswell 1999). By comparison with the line shifts in Table || 
and accepting that several well constrained velocity compo- 
nents will exist in each absorption complex, the level of pre- 
cision achievable with the MM method is a (Aa/a) ~ 10 -5 
per absorption cloud. If systematic effects are to be avoided 
when using a large number (i.e. ~ 50) of absorption sys- 
tems then, from equation Q we require cjq to be known to 
a precision of a(uio) ~ 0.003 cm -1 - an order of magnitude 
more precise than available wavenumber compilations (e.g. 
Morton 1991). 
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Table 2. Shift in the rest frame wavenumber, wavelength and 
velocity space for Act/a = +10 
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AA 
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10 — 2 cm _1 
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AAcr T \9S^ £ i 
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-U.U10 




42 


-u.oo 
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24 
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A 1 tt \1 fi7fl 
I\l 11 AlUlU 




-u. ±o 




A 1 TTT \ 1 S^A 
/\1 111 A-LOO^t 


u.yo 


n '^9 




A 1 TTT \ 1 £fi9 
/\1 111 AlOUZ 


n 

U.^iO 


-U. J.O 


n 09zL 


c; TT \ i croc 

Ol 11 A10ZD 


u. l^t 


fl 039 


n nn^9 
-u.uuoz 


c; tt \i sns 

Ol 11 AlOUO 


1 nfi 

l.UU 


n ^ 
-u.oo 




V_^l 11 AZUOU 


-2 11 


u.oy 


n rj(i 

U. -LOU 


Pv tt \9flfi9 
KjL 11 AZUUZ 


-2 40 


1 02 




Pv tt \9flfifi 
v_ 1 11 AZUUU 


9 ^7 


1.1U 


ft 1 


Tiv> tt X i fins 

re ii Aitjuo 


9 ^7 


n fifi 
-u.uu 


-0 124 


FenA1611 


2.41 


-0.63 


-0.117 


Fe ii A2344 


2.84 


-1.56 


-0.199 


Fe ii A2374 


3.56 


-2.01 


-0.254 


Fe ii A2383 


3.28 


-1.86 


-0.234 


Fe ii A2587 


3.23 


-2.16 


-0.250 


Fe ii A2600 


2.91 


-1.97 


-0.227 


Ni ii A1709 


1.60 


-0.47 


-0.082 


Ni ii A1741 


-1.40 


0.42 


0.073 


Ni ii A1751 


-0.60 


0.18 


0.032 


Zn ii A2026 


4.96 


-2.04 


-0.301 


Zn ii A2062 


3.15 


-1.34 


-0.195 



We therefore require high-resolution, moderate SNR 
QSO spectra and new laboratory measurements of UV res- 
onance transition wavelengths in order to reach a precision 
of o~(Aa/a) ~ 10 -5 for each absorption system. 



5.1 QSO Spectra 

5.1.1 General properties 

All our QSO spectra were obtained at the Keck I 10 m tele- 
scope with the HIRES facility (Vogt et al. 1994). The spectra 
separate conveniently into two samples with several defining 
features. One sample comprises spectra of 28 Mgll/Fell ab- 
sorbers in the spectra of 17 QSOs covering a redshift range 
0.5 < Zabs < 1-8. The second sample comprises 18 damped 
Lyman-a systems (DLAs) in the spectra of 11 QSOs cover- 
ing 1.8 < z a bs < 3.5. This later sample also contains three 
lower redshift Mg n/Fe II systems (i.e. a total of 21 absorbers 
comprises this sample) . Different groups observed each sam- 
ple. Thus, with only a small degree of overlap, we label the 
samples according to the typical redshift of the absorption 
systems. The major difference between the low-z and high- 
z samples is that the low-z sample only contains transitions 
of Mg II and Fe II whereas the high-z sample has a more 
diverse range of transitions present (see Section 5.1.3). 



July. The SNR per pixel ranged from 15-50 with most spec- 
tra having SNR ~ 30. The FWHM was ~6.6 kms _1 (i? = 
45000). The QSOs were generally bright (m v ^ 17.5) and 
so several short (~1000s) exposures were taken of each ob- 
ject. The data were reduced with a combination of iraf[|] 
routines and personal code. Individual frames were overscan 
subtracted, bias frame corrected and flatfielded. Cosmic rays 
were removed by median filtering and the frames were aver- 
aged to form the final image, la error arrays were generated 
assuming Poisson counting statistics. We fitted continua to 
regions of each spectrum containing any of the relevant 
transitions by fitting Legendre polynomials to ~500 kms -1 
sections. Full details of the reduction process are given in 
Churchill (1995, 1997) and Churchill et al. (1999). We have 
also used the spectrum of Q0000— 26 analysed by Lu et al. 
(1996). 

Wavelength calibration of the frames was carried out 
within IRAF: Thorium-Argon (ThAr) lamp exposures were 
taken before and after the QSO exposures and co-added to 
provide a calibration spectrum. ThAr lines were selected and 
centroided to provide a wavelength solution. The solution 
was applied to the corresponding spectra without any re- 
sampling or linearization so that the absorption line shapes 
were not distorted in the process. We also note that no image 
rotator was installed during these observations and so the 
spectrograph slit was not perpendicular to the horizon in 
general. We make a detailed analysis of these points and 
their effect on our results in MOlb. 

The transitions used in this sample are the five Fe 
lines, Fe II A2344-A2600, and the Mg transitions, Mg I A2853, 
MgllA2796 and A2803. The Mg lines have small q r coeffi- 
cients and so act as anchors against which the larger Fell 
shifts can be measured. This set of lines is also very useful 
because it offers a range of line strengths due to the dif- 
fering Fell oscillator strengths (see Table |l]) and, together 
with the similar ionization potentials of Fe II and Mg II, facil- 
itates reliable determination of the velocity structure of the 
absorption complex. In some cases the Mg II lines are satu- 
rated and so the Mgn lines themselves provide only weak 
constraints on Aa/a. However, the weak MglA2853 then 
often serves to provide a good constraint on the velocity 
structure - a further advantage of using Mg/Fe absorption 
systems to probe Aa/a. To illustrate these points, we pro- 
vide an example absorption system, together with our Voigt 
profile fits, in Fig. 1. 



5.1.3 High-z (0.9 < z < 3.5, (z) = 2.1) 

These observations were spread over many observing runs 
from 1994 to 1997. The SNR per pixel ranges from 15- 
40 with most spectra having SNR ~ 30 and FWHM <, 
7.5 kms- 1 (R = 34000). The QSOs themselves are at higher 
redshift and so appear much fainter than those of the low- 
z sample: my ^ 19.0. Consequently, the total integration 
times were much longer than for the low-z sample and many 
more frames were co-added during the reduction process. 



5.1.2 Low-z (0.5 < z < 1.8, (z) = 1.0) 

This sample comprises the same data set anal ysed in W99, 



with some differences as detailed in Section 



6.1 



The ob- 



servations were taken in 1994 July, 1995 January and 1996 



t IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 1. Mg/Fe absorption system towards Q1213— 003 at z = 1.554. The data have been normalized by a fit to the continuum and 
plotted as a histogram. Our Voigt profile fit (solid curve) and the residuals (i.e. [data] — [fit]), normalized to the la errors (horizontal solid 
lines), arc also shown. The tick-marks above the continuum indicate individual velocity components. Note the range of line strengths 
in Fell, facilitating determination of the velocity structure. The large number of Fen transitions and the large number of velocity 
components allows for tight constraints to be placed on Act/ a. In this case, the MglA2853 is also present but is sufficiently weak so as 
not to provide any constraint on Aa/a. In cases where the Mgn lines are more saturated, the Mgl line is typically stronger and can be 
used to constrain the velocity structure. 



An image rotator was installed in 1996 August and so only 
about half the observations were carried out with the slit 
perpendicular to the horizon. 

Most of the data were reduced using makee, the HIRES 
data reduction package written by T. Barlow. This package 
converts the two-dimensional echelle images to fully reduced, 
one-dimensional, wavelength calibrated spectra (calibration 
spectra were taken in the same way as for the low-z sam- 
ple). Some of the spectra were reduced when makee had 
no wavelength calibration facility. In these cases, wavelength 
calibration was carried out using iraf routines. The remain- 
der of the spectra were fully reduced within iraf. la error 
arrays were generated assuming Poisson counting statistics. 
We fitted continua to regions of each spectrum containing 
any of the relevant transitions by fitting Legendre polyno- 
mials to ~ 500 kms -1 sections. Full details of the reduction 
procedures can be found in Prochaska & Wolfe (1996, 1997, 
1999). Outram, Chaffee & Carswell (1999) have also kindly 



provided their spectrum of Q1759+75 which was taken in 
1997 July. 

We also noted a small but important difference between 
the iraf and makee wavelength calibration software, iraf 
makes use of an internal list of vacuum ThAr wavelengths 
and so spectra calibrated with iraf have "correct" wave- 
length scales. However, the makee package calibrates spec- 
tra with a list of air ThAr wavelengths and converts the 
final QSO wavelength scale to vacuum using the Cauchy 
formula for the refractive index (Weast 1979). However, this 
formula does not reproduce the experimental dispersion: it 
shows deviations from experiment at the level of 10 -6 in the 
optical range. We therefore converted the wavelength scales 
of the spectra reduced with makee back to air wavelengths 
using the Cauchy formula and then re-converted to vacuum 
wavelengths with the preferred Edlen (1966) formula. A de- 
tailed discussion of air-vacuum conversion with caveats of 
the above statements is given in MOlb. 
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Figure 2. Heavy element absorption lines from the damped Lyman-a system towards Q2206— 199 at Z — 1.920. We illustrate this 
system since it displays both saturated and unsaturated transitions. However, the metal lines in other DLAs generally have lower column 
densities. The data have been normalized by a fit to the continuum and plotted as a histogram. Our Voigt profile fit (solid curve) and 
the residuals (i.e. [data] — [fit]), normalized to the lcr errors (horizontal solid lines), are also shown. The tick-marks above the continuum 
indicate individual velocity components. In this case, the Fell and Sill lines provide constraints on the velocity structure while the Nin, 
Crn and Znn are strong enough to provide tight constraints on Ao/a. Note that the Cm and Znn A2062 profiles are blended together 
but that the constraints on the velocity structure from all other transitions allow us to make a statistically acceptable fit. 



Theoretically, in damped systems such as those com- 
prising this sample, the combination of lines providing the 
best constraints on Act/ a are those of Znn and Cm. In 
practice, however, these lines are generally quite weak and 
so the parameter constraints are not optimal. Also, the 
Zn II and Cr II A2062 profiles are often blended together (the 
transitions are only separated by 60kms _1 ), reducing the 
strength of the Aa/a constraint. The Ni II lines are also gen- 
erally weak but have the advantage of a variety of widely 
differing qi coefficients. Also note that the only Fell lines 
available for these high-z absorbers are the Fen A1608 and 
A1611 lines. These can be used in conjunction with the two 
Sill lines to provide a reasonable constraint on Aa/a but 
their main advantage is that they provide good constraints 
on the velocity structure. We plot a typical absorption com- 
plex in Fig. 2 to illustrate the above points and also to give 
the reader some idea of the relative strengths of the different 
transitions used in analysing this sample. 



In 6 spectra we have also fitted the two lines of Aim. 
Variations in the incident radiation field could give rise to 
changes in the relative column densities when comparing 
species with ionization potentials above and below the Ly- 
man edge. In our sample, the Al ill velocity structures corre- 
sponded well with those of the singly ionized atoms (consis- 
tent with the conclusions of Wolfe & Prochaska 2000) and 
the value of Aa/a proved to be slightly better constrained 
by including the Aim lines. In some cases, however, some 
velocity components in the Al ill profiles were much weaker 
or stronger relative to the other components when compared 
with the profiles of singly ionized species. In one case, our fit- 
ting algorithm rem ove d many weak Aim components from 
the fit (see Section S.I) and so we removed the Aim data in 
this instance. 
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Table 3. The isotopic structures of the Mg and Si transitions. The 
composite values are given in Table [lj The Mg isotopic spacings 
were measured by Hallstadius (1979) but the spacings for Si are 
estimates based on a scaling of the Mg II A2796 structure. The 
hyperfine structure of 25 Mg is also shown (Drullinger, Wineland 
& Bergquist 1980). The last column shows the relative abundance 
of each isotope (Rosman & Taylor 1998). Isotopic information was 
not used for any heavier ions. 



Transition 


m/amu 


wo / cm 


A /A 


% 


Mtr t X2853 


24 


35051 271 


2852 9636 


79.0 




25 


35051 295 


2852 961 6 


10.0 




26 


35051.318 


2852.9598 


11.0 


Mgn A2803 


24 


35669.286 


2803.5324 


79.0 




25 


35669.300 


2803.5313 


2.7 




25 


35669.310 


2803.5305 


1.5 




25 


35669.360 


2803.5266 


2.8 




25 


35669.370 


2803.5258 


3.0 




26 


35669.388 


2803.5244 


11.0 


Mgn A2796 


24 


35760.835 


2796.3553 


79.0 




25 


35760.853 


2796.3539 


4.2 




25 


35760.913 


2796.3492 


5.8 




26 


35760.937 


2796.3473 


11.0 


SinA1808 


28 


55309.330 


1808.0132 


92.2 




29 


55309.390 


1808.0113 


4.7 




30 


55309.446 


1808.0094 


3.1 


SinA1526 


28 


65500.442 


1526.7073 


92.2 




29 


65500.517 


1526.7055 


4.7 




30 


65500.583 


1526.7040 


3.1 



5.2 Laboratory wavelength measurements 

The laboratory rest wavelengths were measured by various 
groups using Fourier transform spectrometry. We discuss the 
measurements for each ion below, concentrating on those 
issues most relevant to our analysis. We present the wave- 
lengths that we have adopted in Table [i. It is interesting 
to note that, with the exception of the Fell spectrum, the 
wavelengths presented in Table |l| are systematically longer 
than those compiled in Morton (1991) - a widely used list 
of resonance transition wavelengths. For a particular tran- 
sition, the wavelengths are generally consistent but, taken 
altogether, a systematic trend is noted. 

• Mgl and Mgll: These spectra have recently been mea- 
sured by Pickering, Thorne & Webb (1998). The spectra 
of Nil and Nill were excited along with those of Mgl and 
Mg II and several of these lines were used for calibration. The 
Nil spectrum of Litzen, Brault & Thorne (1993) was used 
for reference which itself was calibrated from the Fe I and 
Fen lines of Nave et al. (1991). The MgllA2796 wavelength 
has been measured previously by Drullinger, Wineland & 
Bergquist (1980) and Nagourney & Dehmelt (1981) and the 
wavenumbers are in excellent agreement. Also, a confirma- 
tion of the Pickering et al. (1998) wavenumbers has been 
made and is described briefly in Pickering et al. (2000). 

Since Mg is a relatively light atom, the isotopic structure 
of the line profile must be taken into account. Pickering et 
al. (1998) used the isotopic spacings measured by Hallsta- 
dius (1979) to obtain absolute wavelengths for the isotopes. 
These values are listed in Table |^ and were used in our anal- 
ysis. 



• Sill, Alii and Aim: These spectra have recently been 
measured by Griesmann & Kling (2000). The C IV and Siiv 
spectra were also measured and all spectra were calibrated 
using the Ar II spectrum of Whaling et al. (1995) . At present, 
there are no other similarly precise measurements of these 
spectra. 

To our knowledge, the isotopic spacings for the Si II tran- 
sitions have not been measured. The mass shift and the spe- 
cific mass shift should dominate the volume shift in such 
a light ion. The A1526 transition is very similar to those of 
Mg II and so the specific mass shift should be of the same or- 
der. Thus, to estimate the isotopic spacings in Si II A1526 we 
have scaled the spacings of MgllA2796 by the mass shift: 
AcOi oc LOo/m 2 for Auii the wavenumber shift for a given 
isotope i where m is the atomic mass. We have used these 
spacings to fit a synthetic Gaussian emission line (centred 
at the composite wavelength) in order to estimate the ab- 
solute wavenumbers for the isotopes. We present the results 
in Table ^. We performed similar calculations for the Sill 
A1808 transition. However, the assumption that the specific 
shift for A1808 is similar to that for the Mgll transitions is 
not so justified in this case since the transitions are of quite 
different type (see Table |l|). We discuss any effect this ap- 
proximation may have on Aa/a in MOlb but find that it 
should be negligibly small. 

• Cr II and Zn II: These spectra have been measured inde- 
pendently using two different Fourier transform spectrom- 
eters: one at Imperial College (IC), and another at Lund 
University (LU). The two experiments are described in Pick- 
ering et al. (2000). The wavenumber calibration was based 
on the Fel and II standards of Nave et al. (1991) in both 
cases. 

• Fell: This spectrum, from 1750-3850 A, is presented in 
Nave et al. (1991). The wavelength calibration was done by 
comparison with the Norlen (1973) Aril spectrum. No other 
(similarly precise) catalogues of Fe II lines exist in the litera- 
ture. We have obtained the Fe II A1608 and A1611 wavenum- 
bers from S. Johansson (private communication). 

• Ni II: This spectrum was recently measured by Picker- 
ing et al. (2000). The original spectra were measured dur- 
ing the Mgl and II experiments by Pickering et al. (1998). 
Wavenumber calibration was done using the Ni II standards 
of Litzen, Brault & Thorne (1993) which cover the region 
down to about 2000A. The Nil lines between 1750A and 
2000A were calibrated with the Fel and II lines so as to 
be used as calibration lines for the Ni II spectrum. R. Kling 
(private communication) has confirmed these wavenumbers, 
albeit with lower precision. 

All wavelength calibrations above rely on the the Aril 
spectrum measured by either Norlen (1973) or Whaling 
et al. (1995). However, the Norlen and Whaling calibra- 
tions systematically disagree. The Whaling wavenumbers 
are larger and the difference between them is proportional 
to the wavenumber: 8uj = 7 x 10 -8 u;. Since this difference is 
linear in frequency (velocity) space, as long as we normal- 
ize all measured wavenumbers to the one calibration scale, 
any systematic error will be completely absorbed into the 
redshift parameter when we centroid corresponding velocity 
components in di ffer ent transitions in the absorption sys- 
tems (see Section 3.1 for further explanation). We choose to 
normalize all measured wavenumbers to the Norlen calibra- 
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tion scale and so the wavelengths of the Sill, Alii and Aim 
transitions in Table ^ have been scaled accordingly. 

Also of note is that we only take into account the iso- 
topic structures for the lightest two elements, Mg and Si. In 
the optically thin limit, the centroids of the isotopic struc- 
tures described in Table |^ reduce to the composite values 
given in Table |l| When a line becomes optically thick, dif- 
ferential saturation of the isotopic components causes the 
centroid of the composite line to shift. This effect is dis- 
cussed in detail in MOlb and it is evident that we do not 
need to include isotopic information for the heavier ions. 



6 ANALYSIS AND RESULTS 
6.1 Analysis 

Our analysis technique is based on a simultaneous x 2 min- 
imization analysis of multiple component Voigt profile fits 
to the absorption features in several different transitions. 
Consider a QSO spectrum containing a single velocity com- 
ponent of a specific transition. Three parameters describe 
the Voigt profile fit to such a component: the column den- 
sity, the Doppler width or ^-parameter and the redshift z 
of the absorbing gas cloud. In our case, however, each ab- 
sorption system contains many blended velocity components 
appearing in many different transitions. 

The central assumption in the analysis is that the veloc- 
ity structure seen in one ion corresponds exactly to that seen 
in any other ion. That is, we assume that there is negligible 
proper motion between corresponding velocity components 
of all ionic species. With this assumption comes a reduction 
in the number of free parameters to be varied in a particular 
fit since the redshift of corresponding velocity components 
in different transitions need only be specified once. We dis- 
cuss this assumption and any effect it may have on Aa/a 
in MOlb. 

Further restrictions can be placed on the number of free 
parameters when fitting different transitions simultaneously 
since the 6-parameters of corresponding velocity components 
in different ionic species are physically related. We may write 

°i = -jj- + o turb (6) 

as the 6-parameter of an ionic species i which has a mass 
Mi. The first term describes the thermal broadening within 
a gas cloud which has a kinetic temperature, T, and the 
second term describes an additional turbulent motion which 
affects all ions equally. If we assume that either thermal or 
turbulent broadening dominates for a particular component 
then again we reduce the number of free parameters. If we 
fit transitions of at least two ions with different masses then 
we can determine T and &turb for each velocity component. 

For the present case, we wish to add another free pa- 
rameter to the fit: Aa/a. From equation ^ we can see that 
Aa I a will be constrained by the difference ui z — luo for each 

transition if we fit two or more transitions simultaneously^! 



We have used the program VPFlljJ (Webb 1987) to fit 
absorption profiles to the spectra. Previously, vpfit did not 
allow one to vary T and feturb for each component indepen- 
dently, but dealt only with hi. There was also no provision 
for including Aa/a as a free parameter. We have therefore 
modified vpfit to include T, 6 tur b and Aa/a explicitly as 
free parameters. This differs from the analysis in W99 where 
Aa/a was varied externally to vpfit. The main advantage 
is a significant gain in computational speed, allowing us to 
fit more complicated systems that would have taken pro- 
hibitively long using the method in W99. A minor point 
is that W99 presented 30 absorption systems whereas we 
present 28 using the same QSO spectral data due to slightly 
different subdivision of neighbouring absorption complexes. 

Parameter errors can be calculated from the diago- 
nal terms of the final parameter covariance matrix (Fisher 
1958). The reliability of these errors has been confirmed us- 
ing Monte Carlo simulations of a variety of different combi- 
nations of transitions and velocity structures. 

vpfit imposes a cut-off point in parameter space such 
that very weak velocity components are removed from the fit 
when they no longer significantly affect the value of \ 2 ■ This 
presents a problem for our method since a given component 
may be weak in one transition (and so can be regarded as 
dispensable with respect to y/ 2 ) but may be strong in an- 
other. An example of this can be seen in Fig. 2 where the 
two lowest velocity components do not appear in the weaker 
transitions. Conceivably, dropping even very weak compo- 
nents could affect our determination of Aa/a. Therefore, in 
such cases we observed the trend in the values of (Aa/a)i 
at each iteration i of the minimization routine to see if this 
trend was significantly altered due to line dropping. If com- 
ponents were dropped during a fit then we also re-ran the 
vpfit algorithm, keeping the dropped components by fix- 
ing their column density at the value just before they were 
dropped from the original fit. No cases were found where the 
values of Aa/a from the different runs differed significantly. 

We impose several consistency checks before we accept a 
value of Aa/a. Firstly, the value of \ 2 P er degree of freedom 
must be ~1. Secondly, to check internal consistency, we make 
three different fits to the data with three different conditions 
on the 6-parameters: (i) entirely thermal broadening, (ii) 
entirely turbulent broadening and (iii) a fit where T and 
feturb are determined by goodness of fit. The values of Aa/a 
derived from the three fits should be consistent with each 
other since the choice above should not greatly affect Aa/a 
(i.e. the redshifts of velocity components will not change 
systematically). If this was not the case then the system 
was rejected. Only 2 systems were rejected in this way. Since 
(iii) is the most physically realistic choice we expect y 2 for 
that fit to be less than that for cases (i) and (ii). This was 
indeed the case for all but a few systems. These exceptions 
were found to have very low average temperatures (across 
the velocity structure) or very small values of fc t m-b- The 
difference between y 2 f° r the fit for case (iii) and (ii) or (i) 
respectively was small in these cases. We selected the final 
value of Aa/a from the regime which gave the lowest % 2 
but the results were insensitive to this choice. 



S We require two or more transitions since fitting only one would 
render z and Aa/a degenerate for a single component and nearly 
degenerate for multiple components. 



See 



http://www.ast.cam.ac.uk/~rfc/vpfit.html for details 



about obtaining VPFIT. 
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6.2 Results 

We have analysed the low-z sample in two stages. The first 
was motivated by the fact that we have used different tech- 
niques to those of W99. We took the final velocity structures 
for each absorption system found by W99, used them as first 
guesses for the parameters and determined Aa/a with our 
new methods. In all cases the results were nearly identical to 
those of W99. The second stage was to fit the entire data set 
again, using no information about previous velocity struc- 
tures. These are the results we present below. 

Our results are presented in Table Q We give the ob- 
ject, the QSO emission redshift, z om , the absorption cloud 
redshift, z a bs, Aa/a for that cloud and its associated la er- 
ror. The transitions for each cloud are indicated with the ID 
letters defined in Table |l[ Table ^ shows the weighted mean, 
(Aa/a) w , the unweighted mean, (Aa/a) and the signifi- 
cance level, 8, of the weighted mean for the low and high 
redshift samples and for the sample as a whole. We also in- 
clude the value of the reduced \ 2 > Xrcd (i- e - X 2 P er degree of 
freedom), for each sample where the model is taken to be a 
constant equal to (Aa/a) . 

Our results show that Aa/a deviates from zero at a sig- 
nificance level of 4.1cr over the redshift range 0.5 < z < 3.5. 
We also see a general consistency between the low and 
high-z samples. The values of Xicd m Table |H| are all less 
than unity. This indicates that, if we have miscalculated 
the error bars systematically, we have probably overesti- 
mated rather than underestimated them. This is also re- 
flected in a comparison of the root-mean-square deviation 
from the mean, RMS, and the average size of our error bars, 
(a(Aa/a)): for the sample as a whole, RMS = 1.46 x 1CT 5 
and (a (Aa/a)) = 2.01 x 10 -5 . The RMS is also less than 
(a (Aa/a)) for the low and high-z samples taken alone. As- 
suming we have overestimated the errors on Aa/a then 
we can scale them by a constant factor, S, so as to force 
Xrcd = !• Column 8 in Table |H] shows the resulting increase in 
the significance of (Aa/a) w . Note that the weighted means 
do not differ significantly from the unweighted means for ei- 
ther sample suggesting that the data error bars are reliable. 

We have also conducted an F test on the Aa/a parame- 
ter. We used the final velocity structures from the fitting pro- 
cedure above to re-fit the absorption systems after removing 
the Aa/a parameter. If Aa/a is necessary to describe the 
data then we expect x 2 f° r each system to reduce when we 
include Aa/a. The reduction should be in accordance with 
the statistical significance of Aa/a for each system. We de- 
note \ 2 f° r the fit to system i with and without the Aa/a 
parameter as Xa % an d X%i respectively. If the significance (in 
units of a) of the value of Aa/a for that system is Si, we 
expect that 

M M 

i i 

for M absorbers. 

The 5% can be derived from the last column of Table 
^| Direct calculation of both sides of equation ^yields the 
results in the final two columns of Table [|. ^V(Xai ~X%) <^ 
. 5| for both the low and high-z samples, consistent with 
Aa/a being a required parameter. 

To illustrate the distribution of Aa/a over cosmologi- 
cal time, we plot our results in Fig. 3. The upper panel of 



Table 4. The raw results from the \ 2 minimization procedure. 
For each QSO sight line we identify the QSO emission redshift, 
Zom, the absorption cloud redshift, z a bsi the transitions utilized 
in our analysis and the value of Aa/a for each absorption cloud. 



Object 


2 cm 


Zabs 


Transition" 


Aa/a/10~ 5 






low- 


■z sample 




0000—26 


4.11 


1.434 6 


bcqr 


-0.732 ± 2.283 


0002+05 


1.90 


0.851 


bcnopqr 


—0.239 ± 1.516 


0117+21 


1.49 


0.729 


abccjr 


—0.403 ± 1.298 






1.048 


bcnpr 


— 1.171 ± 1.800 






1.325 


bcpqr 


0.314 ± 0.765 






1.343 


cnpq 


-0.982 ± 0.868 


0420—01 


0.915 


0.633 


abcr 


3.360 ± 5.532 


0450—13 


2.25 


1.175 


bnopr 


—2.590 ± 0.831 






1.230 


bcnpqr 


—1.112 ± 0.798 






1.232 


bcp 


0.611 ± 2.475 


0454+03 


1.34 


0.860 


acnopqr 


0.007 ± 1.266 






1.153 


bcnpqr 


-0.340 ± 2.047 


0823—22 


0.91 c 


0.911 


bcnpqr 


-0.068 ± 0.661 


1148+38 


1.30 


0.553 


bcqr 


-1.882 ± 1.699 


1206+45 


1.16 


0.928 


bcnopqr 


-0.800 ± 1.578 


1213—00 


2.69 


1.320 


abcnopqr 


—1.821 ± 0.893 






1.554 


bcnopqr 


-0.965 ± 0.968 


1222+22 


2.05 


0.668 


bcnpr 


-0.245 + 1.549 


1225+31 


2.22 


1.795 


abcnopr 


—1.473 ± 1.521 


1248+40 


1.03 


0.773 


bcnopqr 


1.398 + 1.260 






0.855 


bcnpqr 


-0.329 ± 1.301 


1254+04 


1.02 


0.519 


abcqr 


-4.486 ± 4.308 






0.934 


bcnpqr 


0.530 + 2.146 


1317+27 


1.01 


0.660 


bcnpqr 


0.348 + 1.781 


1421+33 


1.91 


0.843 


bcnopqr 


-0.111 + 0.793 






0.903 


bcnopqr 


-1.093 + 2.101 






1.173 


bcnpr 


-2.837 + 1.807 


1634+70 


1.34 


0.990 


bcnpqr 


-1.586 + 3.752 






high 


-z sample 




0019—15 


4.53 


3.439 




0.920 ± 3.520 


0149+33 


2.43 


2.140 


dcfghij klnistu 


-3.420 ± 1.800 


0201+37 


2.49 


1.476 


cnoqr 


-0.775 ± 0.878 






1.956 


ehkl 


-0.444 ± 3.243 






2.325 


dcgtil 


3.473 ± 2.414 






2.457 


dgl 


0.902 ± 3.754 






2.462 


ghkltu 


0.520 ± 1.110 


0841+12 


2.20 c 


2.375 


dghijtuvw 


285 + 3 670 






2.476 


dghijklmtu 


-1.351 ± 3.052 


1215+33 


2.61 


1.999 


dcfghij lmtuvw 


-0.972 + 3.411 


1759+75 


3.05 


2.625 


cghklmstu 


-2.199 + 1.521 






2.625 d 


dglmstu 


-0.354 + 1.613 


2206-20 


2.56 


0.948 


bcnpqr 


-2.436 ± 1.724 






1.017 


abcnopqr 


-0.372 + 0.689 






1.920 


dghijklmstuvw 


0.929 + 1.040 


2230+02 


2.15 


1.858 


dghjlnprstu 


-1.615 + 1.699 






1.864 


ghij klmnostuvw 


-1.376 + 0.706 


2231-00 


3.02 


2.066 


ghjlkmtuvw 


-1.590 + 1.265 


2348-14 


2.94 


2.279 


ghl 


-0.244 ± 5.204 


2359-09 


2.31 


2.095 


dghijklstuvw 


-0.946 ± 0.601 






2.154 


dfgl 


-3.101 + 6.211 



"Transitions identified as in Table |lj. 
'This absorber taken from Lu et al. (1996). 

c These QSOs are BL Lac objects and so their emission redshifts 
have large uncertainties. Indeed, the literature values for z cm 
shown here are less than the measured absorption redshifts. 
d This absorber contributed by Outram, Chaffee & Carswell 
(1999). 
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Table 5. Statistics for the two subsamples and the sample as a whole. We give the average redshift, z, for each sample and the number 
of data points, N, contributing to the weighted mean, (Aa/a) , and unweighted mean, (Aa/a) (in units of 10 — 5 ). We also give the 
significance, <5, of the deviation from zero and the reduced \ 2 f° r the weighted mean. S 2 =1 is the significance when we force X^od = ^ 

^rcd 

by scaling all errors on Aa/a by S. The final two columns contain the results of the F test. 



Sample 


z 


N 




(Aa/a) 


S 




5 x 2 =1 

"rod 


S 


V (V 2 - Y 2 1 


y.s 2 


Low z 


1.02 


28 


-0.70 ±0.23 


-0.67 ±0.27 


3.0o- 


0.70 


3.60- 


0.84 


55 


28 


High z 


2.12 


21 


-0.76 ±0.28 


-0.67 ±0.33 


2.80- 


0.68 


3.3<t 


0.83 


37 


21 


Total 


1.49 


49 


-0.72 ±0.18 


-0.67 ±0.21 


4.1(7 


0.68 


4.9cr 


0.83 


92 


49 



IT3 

I 

o 



< 



10 



5 







5 



10 
1 



Redshift 



3 



P- 
1 



"i 1 1 r 



"i 1 1 1 1 1 1 1 1 



H = 68 km/s/Mpc 
(Q M ,n A ) = (0.3,0.7) 
t = 13.86 Gyr 




I 



J I I L 



J I I L 



J I I L 



J I I L 



J I I L 



0.4 0.5 0.6 0.7 0.8 
Fractional look-back time 



0.9 



Figure 3. Aa/a versus fractional look-back time for a currently popular cosmology. The upper panel shows our raw results and la error 
bars: the dots represent the low-z sample and the crosses mark the high redshift sample. Note that the high-z sample does contain some 
lower redshift absorbers. The lower panel shows an arbitrary binning of our results: 7 bins X 7 points per bin = 49 points. The redshifts 
of the points are taken as the mean redshift of clouds within that bin and the value of Aa/a is the weighted mean with its associated 
1<t error bar. 
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Fig. 3 shows the raw values of Aa/a as a function of frac- 
tional look-back time to the absorbing cloud using a flat A 
dominated cosmology (Ho = 68 kms _1 Mpc _1 , Qm ~ 0.3, 
JIa = 0.7). The redshift scale is also given for comparison. 
The lower panel shows an arbitrary binning of the data such 
that all bins have equal number of points (7 bins x 7 points 
per bin = 49 points). We plot the weighted mean for each 
bin with the associated la error bars. At low redshifts we 
see that Aa/a is consistent with zero. Also, the results are 
consistent with a generally monotonic evolution of a with 
redshift. 



6.3 Systematic errors 

The statistical error in our result is now small and our at- 
tention must turn to possible systematic errors. In the com- 
panion paper, MOlb, we explore in detail a broad range 
of potential systematic errors, including the following: lab- 
oratory wavelength errors, wavelength mis-calibration, at- 
mospheric dispersion effects, unidentified interloping tran- 
sitions, isotopic ratio and/or hyperfine structure effects, in- 
trinsic instrumental profile variations, spectrograph temper- 
ature variations, heliocentric velocity corrections, kinematic 
effects and large scale magnetic fields. All but two of these 
are shown to be negligible: atmospheric dispersion effects 
and isotopic abundance variation. In this Section we sum- 
marise these two effects, but first discuss how a general wave- 
length scale distortion affects the low-z and high-z samples 
differently. 



6. 3. 1 Effect of a distortion of the wavelength scale 

To explain our results in terms of an effect other than real 
variation of a, such an effect must be able to mimic a non- 
zero average value of Aa/a for both the high and low-z sam- 
ples alike. Generally, simple systematic effects might cause 
low-order distortions in the wavelength scale (relative to the 
ThAr calibration frames) or may lead to wavelength depen- 
dent mis-centroiding of absorption features. 

The low-z sample (Mgll/Fen systems) is most sensitive 
to systematic effects. The Fen lines have large and positive 
qi coefficients and all lie at lower wavelengths than the Mg n 
anchors. Thus, a systematically negative Aa/a can result 
from a slight 'compression' of the QSO spectra relative to the 
ThAr calibration frames for the low-z sample. However, it is 
more difficult to understand the effect of simple systematic 
errors on the high-z sample for several reasons: 

(i) The transitions involved have a large range of qi coef- 
ficients. In effect, there are three types of transition: those 
with large and positive q\, those with large and negative qi 
and those with small or intermediate values of q\. 

(ii) The different types of transition are mixed in wave- 
length space. In the low-z sample, all the lines with large 
and positive q\ (i.e. Fe II lines) lie to the blue of those with 
small qi (i.e. Mgn lines). This is not the case in the high- 
z sample. For example, two lines with small <ji, Si II A1526 
and Alii A1670, straddle two lines with large, positive qi, 
Fell A1608 and A1611. 

(iii) Different sets of transitions appear in each QSO spec- 
trum. Some spectra may not contain all types of transition. 



The value of (Aa/a) v will be more resistant to system- 
atic errors in the high-z: sample. To illustrate this point we 
have 'compressed' the spectrum of a low-z and a high-z ab- 
sorption system by a factor of 1 x 10 -5 (i.e. the separation 
between any two wavelengths decreases by 1 part in 10 5 ) 
and found a new value for Aa/a. We used the same first- 
guess velocity structures used to derive the results in Table 
|l| and Fig. 3. For the low-z: Mg/Fe system at z = 0.911 we 
find a significant decrease in Aa/a due to the compression: 
Aa/a = (-0.068±0.661) x 10~ 5 -> (-1.106±0.692) x 10 -5 . 
However, for the high-z DLA at z = 2.095 we find a 
small, positive change: Aa/a = (-0.946 ± 0.601) x 10~ 5 -> 
(-0.821 ±0.589) x 10" 5 . 



6.3.2 Atmospheric dispersion effects 

Atmospheric dispersion across the spectral direction of the 
spectrograph slit can be avoided using an image rotator. 
However, as we noted in Sections 5.1. 2| and 5.1.3, many of 
our objects were observed before HIRES was fitted with an 
image rotator. Dispersion across the slit can only lead to 
a 'stretching' of the QSO spectrum relative to the ThAr 
frames. However, at least in the low-z sample, a negative 
value of Aa/a can only be mimicked by a 'compression' of 
the spectrum. The effect of the stretching on Aa/a is not 
quite as clear for the high-z sample (points (i)-(iii) above). 

We quantify the atmospheric dispersion effects in MOlb 
and find that the high-z sample is quite robust against 
this error. However, the low-z results are more sensitive 
to this effect and may have to be reduced from the un- 
corrected value of (Aa/a) w = (-0.70 ± 0.23) x 10" 5 to 
(Aa/a) w ps (-1.52±0.23) x 10" 5 . We stress that the magni- 
tude of the correction here is an upper limit since several ef- 
fects, all difficult to properly quantify, will act to reduce the 
effective compression. Note that since the above correction 
is an upper limit, we do not include any additional errors 
due to uncertainties in our model of atmospheric refraction. 



6.3.3 Isotopic abundance evolution 

We assumed terrestrial isotopic abundances in all species 
when fitting the QSO spectra. However, if the isotopic abun- 
dances in the QSO absorbers are different, small apparent 
shifts in the absorption lines would be introduced, poten- 
tially mimicking a non-zero Aa/a. Aa/a is most sensitive 
to changes in the isotopic abundances of Mgn and Sill 
sinc e these ions have large isotopic separations (see Section 
H]). Timmes, Woosley & Weaver (1995), Timmes & Clay- 
ton (1996) and Gay & Lambert (2000), suggest that the 
abundances of the isotopes 25 ' 26 Mg and 29 > 30 Si decrease with 
decreasing metallicity. The QSO absorption clouds in our 
sample have low metallicity compared to solar (Prochaska 
& Wolfe 1999, 2000; Churchill et al. 1999). Thus, to ob- 
tain an upper limit on the effect this may have had on 
Aa/a, we have removed the 25,26 Mg and 29 - 30 Si isotopes 
from our analysis and re-calculated Aa/a for all the ab- 
sorption systems. Again, we find that the correction is 
towards more negative values, the correc ted va lue being 
Aa/a = (-0.96±0.17) x 10~ 5 (as in Section |6~3~2] we do not 
add any additional errors due to uncertainties in our model 
since the correction we calculate is an upper limit). 
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6.3.4 Tests for unknown, simple systematic errors 

We have also carried out other tests to search for any un- 
known, simple systematic effects. For example, it is unlikely 
that such an effect will be able to mimic the very specific qi 
dependence of the transitions in the high-z sample. Thus, if 
we remove the transitions with large, positive q\ coefficients 
from our fit and find a new value for Aa/a, we expect the 
two values to differ if the line shifts are caused by some 
simple systematic effect. Of course, if we remove transitions 
with large, positive qi, we can only obtain useful constraints 
on Aa/a if anchor lines and transitions with large, negative 
qi are available. Therefore, we have conducted such a test 
on a subset of the high-z sample which contains at least one 
anchor line, at least one transition with large and positive 
qi and at least one with large and negative q\ . We find con- 
sistent weighted mean values of Aa/a both before and after 
line removal. We conducted a similar test by removing lines 
with large, negative q\ and another test by removing the 
anchor lines. We find consistent values before and after line 
removal in each case. Thus, the interpretation that the line 
shifts we observe are caused by a varying a seems robust. 



7 DISCUSSION 

We have conducted an experimental search for variability of 
a over the redshift range 0.5 < z < 3.5 using QSO absorp- 
tion lines as our probe. We used a new method for deriving 
constraints on Aa/a from a comparison of QSO and labora- 
tory spectra which, unlike the previous alkali-doublet (spin- 
orbit splitting) method, exploits the total relativistic effect 
in different atoms (see Section ^|). Comparing the wave- 
lengths of transitions from different multiplets and different 
atoms/ions leads to an order of magnitude gain in precision 
over the alkali-doublet (AD) method. Since the new, many- 
multiplet (MM) method allows us to use all transitions in 
principle, we also gain precision through increased statistics. 

We have analysed the spectra of 49 absorption clouds 
(towards 28 QSOs) and compared the observed and labo- 
ratory wavelengths of those transitions listed in Table |l| 
Our low-z results represent a re-analysis of the W99 data. 
We compare the W99 results with our re-analysis in Table 
and we see that they are consistent. Our results for each 
absorption system are listed in Table ^. Over the entire sam- 
ple we find statistical evidence for time variation of a at the 
4.fcr level: Aa/a = (-0.72 ± 0.18) x 10" 5 (see Table |). 
This result is the 'raw' res ult. Correcting for any possible 
systematic effects (Section 6.3 and MOlb) would increase 
the significance of this deviation of Aa/a from zero. We 
show the distribution of Aa/a with cosmological time and 
redshift in Fig. 3. 

From the analyses we present in this paper and in MOlb, 
we conclude that if our results are due to some effect other 
than a real variation in a, this will only be revealed by fur- 
ther independent observations and analyses using the MM 
method. Optical spectra are very useful for probing the 
range of redshifts we have explored here. However, a sta- 
tistical sample of radio spectra of Hi 21 cm and molecular 
absorption systems has the potential to increase our preci- 
sion limit by another order of magnitude. This effort is ham- 
pered by the fact that very few such systems are known at 



Table 6 . Comparison of results from W99 and the present work 
(low- 2 sample only) in different redshift regimes. Values of Aa/a 
are given in units of 10 — 5 . 







Redshift range 




z < 1.0 


z > 1.0 


0.5 < z < 1.8 


Webb ct al. (1999) 


-0.2 ±0.4 


-1.9 ±0.5 


-1.1 ±0.4 


This paper 


-0.2 ±0.3 


-1.2 ±0.3 


-0.70 ±0.23 



present. A systematic observational search for these systems 
should clearly be made. Also, as the arsenal of 8-10 m class 
optical telescopes with high-resolution spectrographs grows, 
more observations should be devoted to obtaining carefully 
calibrated, high SNR spectra of absorption systems at both 
high and low- z. These will provide a crucial check on our 
results and will allow us to discover or rule out any further 
subtle or unknown systematic effects in the data. 
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